Interference patterns of multifacet 20 x (O-vr-) Josephson junctions with ferromagnetic 
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We have reahzed multifacet Josephson junctions with periodically alternating critical current 
density (MJJs) using superconductor-insulator-ferromagnet-superconductor heterostructures. We 
show that anomalous features of critical current vs. applied magnetic field, observed also for other 
types of MJJs, are caused by a non-uniform flux density (parallel to the barrier) resulting from 
screening currents in the electrodes in the presence of a (parasitic) off-plane field component. 

PACS numbers: 74.50. -|-r, 85.25.Cp, 74.78. Fk 



I. INTRODUCTION 

Multifacet Josephson junctions (MJJs), with critical 
current density jc alternating many times between pos- 
itive j° and negative « — values along the junc- 
tion were intensively treated during the last years. Ini- 
tial studies were motivated by the discovery of high-Tc 
superconductors (cuprates) with d-wave order parameter 
symmetrjsi. In thin film cuprate grain boundary junc- 
tions with a 45° misalignment angle between the two 
electrodes the current jc is changing its sign randomly 
on a scale of a facet. A more controlled MJJ can be pro- 
duced in the form of Nb/cuprate ramp junctions where 
the barrier forms a zigzag line parallel to the cuprate 
crystallographic axes^"-. 

A more general interest in MJJs came from the pos- 
sibility to have a ground state where the Josephson 
phase has a value (p different from both and tt— , 
or to realize tunable plasmonic crystals^^. If a if 
junction is long compared to the Josephson length Aj 
it may carry mobile fractional flux quanta (splintered 
vortices )^^—. A promising option to produce an MJJ 
is given by the Superconductor-Insulator-Ferromagnet- 
Superconductor (SIFS) technologjii^"— , providing expo- 
nentially low damping at low temperatures and a high 
topological flexibility in arranging the and tt seg- 
ments. Recent imaging of the supercurrent distribution 
showed that in SIFSi^ jc is more homogeneous than in 
Nb/cuprate zigzag MJJs^. 

The dependence of the critical current Ic on the ap- 
plied field H for an MJJ is very different from usual 
Josephson junctions. If — — j° and the facets have 
equal size, /c = in zero field but becomes large when 
the field causes constructive interference of the supercur- 
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rents fiowing through the and tt segments. Standard 
calculations of Ic{H) for a junction of length L and width 
W, assuming that the flux density B is homogeneous in 
the tunneling barrier and L,W 4Aj, show that the 
main peaks occur if the flux per facet is $0/2 ($0 is 
the flux quantum), resulting in the maximum current 
Ic — 2LWj^/-K. The other maxima in Ic{H) are much 
lower and symmetric relative to the main Ic maxima^. 

Suppression of /c(0) and appearance of high-fleld main 
maxima were clearly observed for grain boundary MJJs 
in a fleld perpendicular to the substrate plan^^"—. Here, 
randomness in j " , j^ and facet sizes prevent a close com- 
parison to the "ideal" theoretical Ic{H). However, differ- 
ences between experiment and theory also appear for the 
more controllable zigzag MJJs. In particular, the sec- 
ondary maxima between the two main maxima almost 
vanish in experiment and are strongly enhanced outside 
the main peaks^i"— . These features cannot be reproduced 
in calculations by simply taking into account nonunifor- 
mity of the junctions. 

A major step towards understanding Ic{H) was done 
in the framework of nonlocal electrodynamics of MJ Js^-. 
A universal solution has been found for (j){x) in the case 
when the electrodes are formed by an ultrathin film and 
the applied field is perpendicular to them (the definition 
of the coordinates is shown in Fig. [1} . In this geometry, 
the flux density Bz{x) is strongly enhanced compared 
to fioHz, has a maximum in the center of the junction 
and decays to zero towards the edges, see Fig. 2(a) in 
Ref. 22 (the quantity (f>'o{Y) plotted there is proportional 
to the flux density along the junction, i. e. Bz{x) in 
our notation). As a result the main maxima of Ic are 
suppressed compared to the case of uniform B^ while 
the maxima following the main maximum are strongly 
enhanced, see Fig. 4(b) of Ref. [l^. 

In this paper we report the results of our study of Ic (H) 
dependence of rectangular NblAl2O3lNio.6Cuo.4lNb SIFS 
MJJ structures with length L = 200 fj,m and width W = 
10 fj,m. Sections, consisting of 5 /im long O-tt segments 
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FIG. 1: (Color online). Sketch of a SIFS MJJ with 20 facets. 



are repeated = 20 times, as sketched in Fig. [TJ The 
total length L « 3A,/, i.e., these MJJs can be treated as 
short junctions. 
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II. THEORY 

In a short MJJ, see Fig. [TJ the flux generated by tun- 
neling currents is negligible. If By and jc do not depend 
ony, Ic is given by 

rL/2 

Ic — Wmax / {jc{x) sm[(l){x, By) + (j)Q]} dx. (1) 

"^0 J-L/2 

In what follows we consider complex interference pat- 
terns depending on the spatial distribution of flux inside 
the junction. Thus, to show the results of our numeri- 
cal studies, we plot Ic/Ico as a function of H/Hq, where 
Ico = Jc^W and = ^q/iiqAL. A is the (average) 
effective magnetic junction thickness accounting for flux 
penetration into the superconducting electrodes and the 
effect of screening currents. (A w 2Al for a junction in a 
bulk sample, where is the London penetration depth); 
below we will also allow A to be different in the and tt 
parts. 

The phase (f>{x, By) is obtained from the equation 



d0 
dx 



2tt By {x)K 



(2) 



The calculated Ic{H) for constant A is shown in Fig.[2][a). 
The two main maxima of Ic occur at \H\/Hq = N. The 
next maxima are located at \H\/Hq — 3N. Further, 
4 = at \H\/Ho = 2N. 



FIG. 2: (Color online), (a) Calculated critical current of a 
20 X (O-TT-) MJJ as a function of applied field, for j° = 
facets of equal size and a homogeneous field B || y. (b) A 
"typical" experimental pattern (T = 4.2 K), when H is mis- 
aligned relative to the substrate plane {9 ~ — 1°)- 
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FIG. 3: Ic{H) curves for (a,b) parallel {6 = 0°) and (c,d) 
perpendicular (9 = 90°) field orientation. Inset in (c) is a 
sketch of B in the (y, z) plane for 6 = 90°. Inset in (d) shows 
the profile By{x) used to calculate Ic{H) for 9 = 90°. 



III. EXPERIMENT 

The current-voltage curve of the MJJs that we stud- 
ied experimentally (T k, 4.2 K) are non-hysteretic. We 
determined Ic by using a voltage criterion Vc-c — 1 /iV ex- 
cept for the data shown in Fig. [2][b) , where we have used 
a different measurement setup, with T4r = 0.2 /iV. Thus, 
we have an Ic detection limit (parasitic Ic background) 
/„ = V„/R (R ~ 0.05 rj). /„. w 4^A in Fig. ^h) and 



w 20 /iA in all other experiments. When the sample was 
aligned "by eye" to have H \\ y, the IdH) patterns (see, 
e.g., Fig. ^h); in comparison to Fig. ^a) the applied 
field scale roughly corresponds to ±30i?o) resembled the 
ones of zigzag junctions^i^ (where H \\ z), but not the 
calculated "ideal case" pattern shown in Fig. [2][a). In 
particular, the measured IciH) shows a small local max- 
imum at _ff = and high maxima following the main 
maxima. 

The model considered in Ref. [l^ (ultrathin film grain 



3 



boundary MJJ in the (x, y) plane) cannot be applied, 
since (a) the geometry of SIFS MJJ is more complicated 
and (b) the validity conditions for the nonlocal model are 
not satisfied. However, for the SIFS (and also zigzag) 
MJJs an applied field ff || i should be subject to field 
focusing effects, with a maximum of By{x) in the center 
of the junction^^i^i. In our MJJ aligned "by eye" a com- 
ponent Hz may arise from a slight out-of-plane misalign- 
ment of the applied field. This small field leads to strong 
modifications of Ic{H). To test this idea we measured 
Ic{H) for different angles between H and the substrate 
plane (see Fig. [!}, using two perpendicular coils oper- 
ated in linear combination and creating Hy = H cos 9 
and Hz = Hsm9. 

Fig. [31(a) shows Ic{H) measured at an angle, which, 
after having measured and analyzed the angle dependent 
Ic{H) (see below), has been identified to be very near 
6* = 0°. The shape of this curve is much closer to the 
"ideal case", although differences occur. First, there is 
an /c maximum aX H = 0. This naturally arises, because, 
typically, | | ^ and thus the cancelation of the super- 
currents carried by the and tt parts is incomplete. Sec- 
ond, an Ic peak appears at H/ Hq = 40 = 2N. If both the 
and TT segments had the same length and the same A, 
at this field Ic should cancel within each segment. How- 
ever, A may slightly differ for the and tt parts. Then, 
Ic cannot cancel at the same field in both the and the tt 
segments and a peak may appear in Ic{H). Further note 
that the measured Ic{H) is asymmetric with respect to 
positive and negative H . This effect arises from asym- 
metries in the ferromagnet layer magnetization which we 
do not take into account to keep the discussion short. 

Taking asymmetries of jc and different values Ap and 
A^ of the and tt parts into account, using Eq. ([T]) 
and Eq. we calculate the Ic{H) pattern shown in 
Fig. ISKb). In Fig. [3Kb) we have used A^/Aq = 1.38 
and jc/jc — —0.8 to achieve the best fit to the exper- 
imental curve. This translates to j° « 42A/cm^ and 
w — 34A/cm^ in physical units. These parameters 
are reasonable and the resulting Ic{H) describes the data 
well. Also note that the ratios Jc/jc and A^/Aq can be 
determined independently and in a straightforward way, 
by analyzing the Ic maxima at, respectively, H — and 
H — 20Hq. Comparing the field axes of Figs. [31(a) and 
(b) we obtain Aq « 134 nm and A^r ~ 186 nm in reason- 
able agreement with the value of 2Xl{Xl ~ 90 nm for 
Nb). Further, using 



measurement 



calculation 

8 = 6.21 (b)^ 



X'j = y/<i>o/2TTHo\jl\d', with i = 0, TT. 



(3) 



Here fiod' is the inductance per square (with respect to 
in-plane currents) of the electrodes and d' w A for elec- 
trode thicknesses Al, we find a normalized junction 
length I = L/2X^j + L/2A} w 3.2, i.e., this MJJ is in the 
short limit. 

Having analyzed Ic{Hy) and fixed the junction param- 
eters j7> ■'^o and A^ we turn to the other limit of 
H II z, {9 = 90°). The corresponding measured Ic{H) 
pattern is shown in Fig. [3Kc). Its shape resembles the 




FIG. 4: Ic{H) patterns for different misalignment angles 9. 



ones of zigzag junctions^ — or of the SIFS MJJ with H \\ y 
aligned "by eye", c.f. Fig. [21(b). Note that our MJJ is 
sensitive to Hz, because screening currents flowing in- 
side the electrodes create a non-uniform field component 
By[x) inside the junction barrier, as sketched in the in- 
set of Fig. [SJc). To calculate Ic{Hz) we need to find 
a profile By{x) caused by the applied field H \\ z and 
solve Eqs. (jlj and (0) using this By{x). Similar to the 
case of a junction having overlap geometry^^i^i and to 
the nonlocal grain boundary junctioi>2^ we expect By (x) 
to have a maximum in the center of the MJJ (x — 0). 
It can be approximated bl^^^— By{x) oc cos{ttx/L) if 
the junction is "naked", i.e., without idle regions sur- 
rounding the Josephson barrier. In our case there are 
5 fim wide idle regions. Then, By will be nonzero at 
X = ±L/2. In addition, asymmetries in By{x) may oc- 
cur due to asymmetries in the junction layout. Thus, we 
approximate By{x) by a Taylor expansion near x = 0, 
i.e., By{x) = —f^oHz{l + oi^ -I- 02^^), where ^ — 2x/L 
and -1 < ^ < 1 . Solving Eq. we get the Ic{H) 
pattern shown in Fig. [3Kd). We have used ai — —0.1 
and 02 = —0.45 to reproduce the data. While the over- 
all shape of Ic{H) strongly depends on the value of 02, 
the nonzero choice of ai is somewhat cosmetic and was 
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merely necessary to suppress a "shoulder" in the calcu- 
lated Ic{H) oscillations appearing for flux values above 
the main /c maximum. For lower values of H this pa- 
rameter has only a small effect on Ic- The resulting field 
profile is shown in the inset of Fig.[3][d). Further, by com- 
paring the abscissas of the calculated and experimental 
curves, we find a field focusing factor / « 23, which is a 
reasonable value for the large MJJ we study. 

Having reproduced Ic{H) for both H \\ y and H \\ z, 
we consider the general case where H is tilted by an angle 
9. No additional parameters are required for calculating 
Ic{H). We obtain 

ByiO = -^loHiif sme-cose) + {aiC+a2^^)f sine]. (4) 

Similar to the geometries analyzed in Refs. [2^1251 at a 
"dead angle" 9d — arctan(l//) the uniform field term 
vanishes and one obtains an anomalous IdH) depen- 
dence with diverging period. The linear and quadratic 
terms in By{x) can have different signs relative to the 
constant term depending on 9. For angles between 9^ 
and 0°, and for a negative coefficient 02, By(x) peaks at 
the junction edges while at all other angles the maximum 
in By is in the junction center. 

In Fig. |4]we show a sequence of Ic{H) patterns start- 
ing from 9 = Q.2 to 9 = —2.43° comparing the mea- 
sured Ic{H) with calculated curves. All data have been 
taken in a single cooldown. In the simulated curves the 
junction parameters and coefficients 01,02 are the same 
as for the cases 9 = {) and 9 = 90° discussed above. 
The patterns taken a.t 9 = 6.2°, see Fig. Hl^a) differ only 
marginally from the case oi 9 = 90°, shown in Fig. ^a). 
The main difference is the field scale due to the s'm9- 
reduced By[x) field caused by H^. In particular, for the 
cases 9 = 6.2° and 9 = 4.5° the enhancement of the 
higher order Ic maxima (relative to the main maximum) 
can nicely be seen. When approaching 9d{~ 2.5° for our 
sample), Ic{H) stretches anomalously with an increased 
amplitude of the low-field Ic oscillation and with an in- 
creased field modulation period. Further, for angles be- 
tween 9d and 0° (where, according to Eq. (j?]), the field 
profile By{x) reverses shape) the lower order Ic maxima 
become enhanced while the higher order maxima are sup- 
pressed. These features are well reproduced by the calcu- 
lated curves IdH). Ic{H) at = 0° is shown in Fig.[31[a). 
For negative values of 9 the shape of Ic{H) rapidly devel- 
ops into the pattern found for || z, without anomalies 
at 6* = -9d. 

A "dead angle" has been found for junctions of var- 
ious geometries^^J^. To further test this for SIFS junc- 
tions having large idle regions we have measured Ic{H) 
of a 50 /im x 10 /im coupled SIFS junction, for H ap- 
plied under different angles 9. Fig. [5] shows the resulting 
curves for three values of 9. For = 0° and 90° Ic{H) 
is close to a Fraunhofer pattern, with a compressed field 
scale for the 9 = 90° case; by comparing the modulation 
periods of these two curves we find / « 2.3 and calcu- 
late 9d ~ 23.5°. The dependence IdH), measured at 
9 = 23.65° and shown in Fig. [51 indeed is by no means 
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FIG. 5: (Color online). Ic(H) curves for a 50 x 10^x0^ SIFS 
junction for 9 = 0°, 90° and 23.65° (^ Od). 



Fraunhofer like. 

The value of 9d depends on the lead geometry via the 
field focusing factor / which is large for long junctions 
with wide leads. While 9d « 23.65°, found for our 
junction, is not very critical if the sample is aligned "by 
eye" , other geometries, particularly in the context of long 
junction physics may have larger values of / and corre- 
spondingly lower values of 9d. 



IV. SUMMARY 

We have studied a 20 x (O-tt-) multifacet SIFS Joseph- 
son junction. Its IdH) dependence, when the applied 
magnetic field is aligned accurately parallel to the junc- 
tion plane, can be described well by the standard lin- 
ear phase ansatz (constant field), taking into account 
lie I 7^ be I ^-iid different values for the effective mag- 
netic junction thickness A in the 0- and 7r-parts. On 
the other hand, a variety of IdH) patterns can be ob- 
tained as a result of a small misalignment and field focus- 
ing. When the perpendicular field component dominates, 
the IdH) patterns are similar to the ones measured for 
Nb/cuprate zigzag O-tt junctions. Specific features are: 

(a) a suppressed amplitude of the main Ic maxima and 

(b) a strong asymmetry between the lower order and 
higher order maxima. If the field is off-plane one may 
meet a "dead angle" 9d where the flux density By caused 
by Hy is cancelled by the constant part of a non-uniform 
By{x) induced by Hz. The effect is not restricted to mul- 
tifacet SIFS junctions. If the sample is aligned so that 
9 K 9d, Id-^) is anomalous and may lead to erroneous 
conclusions about the sample quality or even the physics 
investigated. 

Although the MJJ investigated here has a too large 
jc asymmetry to from a tp junction, the effects discussed 
here should be taken into account for the right interpre- 
tation of experimental results obtained for such devices 
in applied magnetic field. 
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